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I. INTRODUCTION 


One of classical problems in astronomy and physics is the three-body problem in Newto¬ 
nian gravity (e.g., |lH3[|). The gravitational three-body problem is not integrable by analyti¬ 
cal methods. As particular solutions, however, Enler and Lagrange found a collincar solution 
and an equilateral triangular one, respectively. The solutions to the restricted three-body 
problem, where one of the three bodies is a test particle, are known as Lagrangian points 


L \, L/ 2 i L 3 , L 4 and L 5 


a 


Lagrange’s equilateral triangular solution has also a practical 
importance, since L 4 and L 5 for the Sun-Jupiter system are stable and indeed the Trojan 
asteroids are located there. Even for the Sun-Earth system, asteroids were also found around 

n 

L 4 by recent observations [4J. 

Recently, Lagrangian points have attracted renewed interests for relativistic astrophysics 
5fjip!, where they have discussed the post-Newtonian (PN) corrections for Lagrangian points 
6 ], and t he g ravitational radiation reaction on L 4 and L 5 analytically [f] and by numerical 


methods 8 -jlO]. It is currently important to reexamine Lagrangian points in the framework 
of general relativity. As a pioneering work jll], Nordtvedt has pointed out that the location 
of the triangular points is very sensitive to the ratio of the gravitational mass to the inertial 
one, though his analysis does not fully take account of the first post-Newtonian (1PN) 
terms. Along this course, it might be important as a gravity experiment to discuss the 
three-body coupling terms in the PN force, because some of the terms are proportional to 
a product of three masses such as m 4 x m 2 x m 3 . Such a triple product appears only for 
relativistic three (or more) body systems but it does not for relativistic compact binaries 
nor Newtonian three-body systems. In addition, it has been pointed out that three-body 
interactions might play important roles for compact binary mergers in hierarchical triple 
systems 12 16]. Very recently, moreover, a first relativistic hierarchical triple system has 


been discovered by Ransom and his collaborators 




For three finite masses, in the 1PN approximation, the existence and uniqueness of a 


col 


inear solution corresponding to Euler’s one have been shown by Yamada and Asada 


18l. Il9| . Also, Ichita et ah, including two of the present authors, have shown that an 
equilateral triangular solution is possible at the 1PN order, if and only if all the three 


masses are equal 


20]. Generalizing this earlier work, Yamada and Asada have found a PN 


triangular equilibrium solution for general masses with 1PN corrections to each side length 
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21|. This PN triangular configuration for general masses is not always equilateral and it 
recovers the previous results by Krefetz [5|] and Maindl 6 ] for the restricted three-body case. 

In Newtonian gravity, Gascheau proved that Lagrange’s equilateral triangular configura¬ 
tion for circular motion is stable [ 22 |, if 


mim 2 + m 2 m3 + m^mi 1 
(mi + m 2 + m 3) 2 27 


( 1 ) 


Routh extended t 
for stability as [23] 


re result to a general law of gravitation oc 1 /r k , and found the condition 


mim 2 + m 2 m 3 + m^rrii 1/3 — k 
(m 1 + m 2 + m 3 ) 2 3 \ 1 + k 


( 2 ) 


For the restricted three-body limit as m 3 —> 0, Douskos and Perdios examined the stability 
of L 4 (L 5 ) in the 1PN approximation, and they obtained a region of stability as 24] 

m 2 17\/69 

-i-< l*o - 7^~ e ’ 3 

mi + m 2 486 

where we assume mi > m 2 without loss of generality, the Newtonian value fio = (9—\/69)/18, 
and we define 

GM 


£ 


r'-'i 


( 4 ) 


for the total mass M and each side length t of the Newtonian equilateral triangle. Singh and 
Bello also discussed the stability in the restricted three-body problem by taking account of 
not only the general relativistic effects but also radiation pressure and a small perturbation 


in the centrifugal force 


25 


26]. For three finite masses, however, it is not clear whether 


the PN triangular configuration is stable. In this paper, we study this issue by fully taking 
account of all the 1PN terms. 

This paper is organized as follows. In Sec. im we briefly summarize the PN triangular 
equilibrium solution for three finite masses. In Sec. IIII1 we consider the stability of the 
solution at the 1PN order. Section m is devoted to the conclusion. Hereafter, we take the 
units of G — c — 1. 


II. A POST-NEWTONIAN TRIANGULAR SOLUTION 


In this section, following Ref. 2l|, we summarize a derivation of the PN triangular 
solution to the relativistic three-body problem for general masses. We take account of the 
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terms at the 1PN order by employing t 


in the standard PN coordinates as 


27 


re Einstein-Infeld-Hoffman (EIH) equations of motion 


30] 


d 2 r 


m K - 


K 


dt 2 


£ 


r A K- 


m K m A 


A^K 

,2 


' AK 


4 E 

B^K 


m B 


EErl 1 

C^A 


rAK ■ r CA 

9 r 2 
zr CA 


3 

+v 2 k + 2v\ - Av a ■ v K - - (v A ■ n AK f 


~ ( Va ~ Vk ) 

A^K 


m K m A 


1 AK 


/ q . \ | 7 \ ^ \ ^ m K rn A rn c , . 

n AK ■ (3 v a - 4n A 0 + ^ r ° A r -’ ( 5 ) 

'AKiCA 


A+K C^A 

where rj and Vj denote the location and the velocity of each body in an inertial frame and 
we define 


ru = r T - rj, 


ru = \ru |, 
Tu 


nu = 


ru 


( 6 ) 

( 7 ) 

( 8 ) 


In the following, we assume circular motion. 

We consider a PN triangular configuration with 1PN corrections to each side length of a 
Newtonian equilateral triangle, so that the distances between the bodies are 


r LJ = e(l + pu), 


(9) 


where /, J = 1,2,3 and pu{— pji ) is dimensionless PN corrections (see Fig. [TJ) . Because 
of the circular motion, £ and pu are constants. Note that we neglect the terms of second 
(and higher) order in e henceforth. Here, if all the three corrections are equal (i.e. p ±2 = 
P 23 = P 31 = p), a PN configuration is still an equilateral triangle, though each side length is 
changed by a scale transformation as £ —» £{1 + p). Namely, one of the degrees of freedom 
for the PN corrections corresponds to a scale transformation, and this is unimportant. In 
order to eliminate this degree of freedom, we impose a constraint condition 

r'12 + r 23 + r 31 
3 

which means that the arithmetical mean of the three distances of the bodies is not changed 
by the PN corrections. Namely, 



P12 + P23 + P31 — 0 . 


( 11 ) 
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Please see also Ref. 21] for imposing this constraint. 

Let uji denote the angular velocities of the Jth body with PN corrections. The EIH 
equation of motion for m\ becomes 

1 


2 M _ 3 M 

- u i r ni = --7y r Ni + o E ihi£ - bTT r-o— 

P 2 t ^/vl + z/ 2 z / 3 + v- 


x {[v 2 {v\ —V 2 - l)pl2 + ^ 3(^1 - "3 — l)p3l]ni + VSu 2 U 3 (p 12 - P3i)n_Li}, (12) 

where the mass ratio vj = m//M, t*n/ and v^i are the Newtonian location and velocity, 
respectively, rij = r N //r N /, and n LI = « N //% i- Note that rij _L n LI in the circular motion. 
The PN term <5 eihi is defined by 

. 1 M 1 

Orth 1 — 


16 P a/z/I + v 2 v 3 + z/f 

x / {16 (z^2 + ^ 2^3 + ^ 3 ) [3 — (^ 1^2 + ^ 2^3 + ^ 3 ^ 1 )] 

+ Qv 2 is 3 [2(z/ 2 + v 3) + z^2 + 3 u 2 u 3 + ^3 ]} ^ 1 

+ 3\/3u 2 u 3 (u 2 « Z2 3 )(5 - 3z/ 1 )n_Li|. (13) 

One can obtain the equations of motion for m 2 and m 3 by cyclic manipulations as 1 — > 2 — y 

3 -> 1 

The PN triangular configuration becomes an equilibrium solution in the circular mo¬ 
tion if and only if the following conditions (a) and (b) simultaneously hold: (a) the term 
proportional to n±j vanishes and (b) all the angular velocities are the same. 


As a result, we obtain the PN corrections that satisfy the above conditions as 21] 


Pi 2 = 24 [(^2 - Rs)(5 - 3zq) - (z / 3 - i/i)(5 - 3z/ 2 )]e, 

P 23 = ^[(p? - ^i)(5 - 3 z/ 2 ) - (zq - u 2 ){5 - 3v 3 )]e, 

P 31 = ^[(^1 - P 2 )(5 - 3z/ 3 ) - (z / 2 - z/ 3 )(5 - 3z/i)]e, 


(14) 

(15) 

(16) 


which give a PN triangular equilibrium solution for general masses. In this case, the angular 
velocities of the bodies are equal (i.e. oq = oj 2 = co 3 = lo). These corrections for the 
restricted three-body problem as v 3 —> 0 reduce to previous results Ik- 

Substituting Eq. (jT4|) and Eq. (1T6|) into Eq. (1T2|) . we obtain the angular velocity of the 
three bodies 


U — CUn(1 + CUpn), 


(17) 
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where the Newtonian angular velocity cun 


yjM /£ 3 and the dimensionless 1PN correction 


^PN = -77(29 - 14(z/!Z/ 2 + z/ 2 z/ 3 + i/ 3 i/i)]e. 
lb 


(18) 


Using a relation zq + + zz 3 = 1, we can show cDpn < 0, that is, u < cun for the hxed system 

parameters t and M. In other words, the PN triangular configuration is always smaller than 
the Newtonian one if the masses and angular velocity of the three bodies are hxed. 


III. THE STABILITY OF THE POST-NEWTONIAN TRIANGULAR SOLUTION 


Next, we study the stability of the PN triangular solution by taking account of linear 
perturbations in the orbital plane. It is convenient to use the corotating coordinates with 
the origin as the center of mass even after adding perturbations. Therefore, the number 
of the degrees of freedom for the perturbations decreases from six to four. One of them 
corresponds to a perturbation in the angular velocity, and the three other perturbations 
denote changes in the shape and size of the PN triangle. 

We consider two perturbations corresponding to changes in the distances r \2 and r 3 1 as 

ri2=i(l + P12 + SX12), (19) 

r 31 =£(l+p 31 +5x3i), (20) 

where X 12 and X 31 are perturbations in the distances and 5 is a bookkeeping parameter that 
denotes the smallness of the perturbations. And, a perturbation ^23 denotes a change in 
the angle tp 23 between r 32 and r 33 as 

<^23 = ^ + \/3)023 + ^23- (21) 

These three perturbations mean changes in the shape and size of the PN triangle. For the 
remaining one of the degrees of freedom corresponding to a change in the angular velocity, 
we denote 


$12 — ©12 + b ( J , 


( 22 ) 


where 612 and a denote the direction of 772 to the reference frame and a perturbation in it, 
respectively. @12 is the unperturbed direction, which satisfies the equation as 

d@i2 


dt 


— <^n(1 + ^pn). 


(23) 
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Note that the 1PN corrections pu satisfy Eq. (fill) . while the perturbations are arbitrary. 
Figure [2] shows a schematic figure for these four perturbations. The above choice of the 


pertur 


27 


29 


Dations is convenient, because we can avoid directly using the PN center of mass 


Note that the perturbations X 12 , X 31 , ^23 , and a have not only the Newtonian terms but 
also the 1PN ones. For instance, the perturbation cr can be expanded as 


a — (Tn + crpN, 


(24) 


where <7n and <Tpn(= 0(e)) are the Newtonian term and the 1PN one, respectively. In the 
following, we neglect the terms of second (and higher) order in <5. Namely, we calculate to 
the terms of order e x 6 (i.e. the linear perturbation at the 1PN order). 

Defining a new variable as X = y 31 — y 12 , we obtain the equation of motion for r 12 (see 
Appendix [A] for a detailed derivation of the equations). Its radial part is 


9 3\/3 

(D 2 - 3) X i2 ~ 2 Da - -u 3 X - -^-^ 23 


+ s 


— <1 4a/3(^! - z/ 2 )(7- 9v 3 )v 3 D 


+ (36z/| + 234zqz/ 2 - 146^2 + 261 u 2 u 2 - 488z/it/ 2 + 155u 2 + 63z/ 3 - 155z/ 2 + 137 v x 
- 585)|xi 2 - ^(27z/ 3 + 135zqz/ 2 - 21u 2 + 135u 2 z/ 2 - 210z/in 2 + 24z/ 2 + 27z/ 3 - 21^ 

+ 24iq — 155 )Da - 7/3 j 4v / 3(9iqi/ 2 + 10z/ 2 + 9u 2 — Qui — 4 )D — (216z/ 2 + 288zqz/ 2 

32 ! 


154t/ 2 + 171z/ 2 — 38z/i + 420) >A" H-z/ 3 <( 4(18 z/ 2 + 2 l 7viv 2 — 2v 2 + 9v 2 + 14zq 


J 32 

- 12) D + v^3(51 u 2 + 1Uuiu 2 + 2u 2 + 87n 2 - 120z/i + 155) U 23 


= 0, 


(25) 
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and the tangential part is 


3^ 


9 


2 £>Xi 2 + D~a - -v-iX + 7^23 


+ £ 


32 


4(9^ + 45z/iz/2 + 92/ 2 + 452/ 2 2/ 2 


— 30z/iZ/ 2 — 18i/ 2 + 9Z/J 5 + — 18 2/1 + 61 )D + 3^2/3(122/3 — 6//1 // 2 + 14z/ 2 — 152/ 2 

+ 4 2/1 — 5) l’Xi2 — — “f (32/| + 122 /i 2/ 2 — 182/ 2 + 32/j — I82/1 + 10)14 2 — 3 v / 3(2/i — 2/ 2 ) 


X 2/3(92/2 + 92/! + 4)14 


} a + 


4(182/3 + 272/! 2/ 2 + 82/ 2 + 92/ 2 + 16 2/1 — 12)14 


V / 3(362/ 2 + 722/12/2 — 542/ 2 + 8l2/j" — 9O2/1 + 160) ^ A + — 2/ 3 -^4v/3(92/i2/ 2 + 


82/2 


+ 92/ 2 — 4)14 — 9(2l2/ 2 + 142 /i2/ 2 — 102/ 2 + 132/ 2 — 82/1 + 45) >/2 23 


= 0, 


(26) 


where 14 denotes a differential operator with respect to a normalized time 1 = 

In the same way, we obtain the equation of motion for 7*31 and its radial and tangential 
parts are 

/ „ 9 \ ( 3, a/3 \ 1 

+ £ 


! - 3 + jU 2 'J X - (2 D + ~~~2/ 2 | l \) 2 3 


--4^3 
32 ’ 


(Z) 2 — 3)xi2 — 2Z)cr + ^.D 

x (2/3 — 2/1)(7 — 92/ 2 )2/ 2 14 + (362/3 + 23427]2y 2 — 146273 + 2612/ 2 27 3 — 4882/12/3 + 15527 3 
+ 6327^ - 15527 2 + 137271 - 585) jxi2 - ^-(272/f + 13527^ - 21i/f + 1352/ 2 2/ 3 - 2IO2/12/3 

+ 24273 + 27 uf - 2127 2 + 24271 - 155 ) 14(7 - — j 4V327 2 (927 3 2 + 9271Z7 3 + 8273 - 4271 - 4)14 

32 I 

— (I8O273 + 270271 Z7g — 22427^ + 19827 2 Z7 3 + 82/1273 + 4192/3 + IO82/1 — 542/ 2 + 3212/1 
+ 165) jx + ^|4(27i/| - 392/f - 272/ 2 2/ 3 + 1652/! 2/3 - 542/ 3 + 362/ 2 - 1022/1 + 191)14 

A 1 


3^32/2(512/f + II42/12/3 + 22/3 + 872/ 2 - 1202/1 + 155) 1^23 


= 0, 


(27) 



2D X i2 + D 2 a+[2D- ^-u 2 j X + (d 2 - ^ 2 ) V 23 


+ £ 


32 


4 ( 9^3 + Abvivl 


+ 91/3 + 45^1/3 — 30z/ii/ 3 — 18 z/ 3 + 9Z/J 5 + 9z/ 2 — 18z/i + 61)£> — 3V3z/ 2 (12£/ 2 — 6//1 


^3 


+ 14z/ 3 - 151/3 + 4z/i - 5) |xi 2 - — j(3z/f + 12*4 ^3 - 18i/ 3 + 31/3 - I 8 Z /3 + 10)L > 2 

- 3 v / 3 (z /3 - z/i)(13 - 9l , 2 s )v 2 D |cr + ^| 4 ( 9z/ I “ 19t/ 3 _ 9i/ 2 z / 3 + 27i/iz/ 3 - 2z / 3 - 2z / 2 

- 10 //, - 49)L> + V3(72 v\ + 54z/iZ / 3 - 12 z / 3 + 36//( - 78^ + 1A5)u 2 \x - ^|4(3z/| 
+ 12 z/ 3 i / 3 — 18z/ 3 + 3i / 2 — I 8 Z /3 + 10).D 2 — 12\/3(9i/| + 9 z/iz / 3 + 12z / 3 — Au\ — A)u 2 D 


— 27(21z/ 3 + I 41 / 3 I /3 — lOz /3 + I 3 Z /3 — 8 zq + 45)z/ 2 z'023 


= 0 . 


(28) 


First, we study the condition for stability in Newtonian gravity. In the Newtonian limit 
e —> 0 , the equations of motion (f25l) - (|28l) for the perturbations are rearranged as 


/ 


(L> 2 -3) -2D 
2D D 2 


9 

r 

3\/3 


~r 3 


3^3 


- 1/3 

(L> 2 -3) -2D D 2 -3 + -v 2 


\ 


2D 


D 2 2D- 


3V3 


-v 2 


-Vi 

u 
9 

~A 3 _ 
2D + ^i / 2 

H 2 - ^2 


\ 




/ 


(7 




^23 ) 


= 0 . 


(29) 


Carrying out a classical stability analysis (e.g. Q), we obtain the eigenvalue equation, so 
called the secular equation as 


, 2,32 ,,3/3 2 , 1 +V1-27C\ / 2 , 1 - VI - 27C 


A (A + 1) A + 


A + 


= 0 , 


(30) 


2 ) \ 2 
where A is the eigenvalue and V = viv 2 + u 2 u 3 + z/ 3 z/ 3 . Since Eq. (130]) has only the terms of 
the even-order of A, if Re(A) V 0, both the positive and negative real parts are allowed and 
then the equilibrium configuration is unstable. In fact, the roots of Eq. (]30l) are 


\ n 3 , • 3 , , 1 + VI - 27E . , / 1 - VI - 27C , Q1 . 

Ao = 0, Ai ± = ±i, \ 2 ± = ±\l ---, A 3± = ±y---• (31) 

Therefore, Lagrange’s equilateral triangular solution is stable if and only if A 2 ± and A 3 -t are 
purely imaginary. Hence, it is necessary and sufficient that 


1 - 27E > 0. 


(32) 
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This is nothing but the Newtonian condition Eq. ([I]) for stability of Lagrange’s solution. 

Next, let us consider at the 1PN order. The EIH equations of motion (t25|) - (1281) for the 
perturbations are 


where 


^M\ 1 A/13 Mi 4 ^ 



CO 

mn 

CN 

MC 


a 

M31 m 32 m 33 m 34 


X 

\M 4 i M 42 M43 M 44 ) 


\'fe/ 


= 0, (33) 


Mn = (77 2 — 3) — — £ |4\/3(zq — z/ 2 ) (7 — 9u 3 )u 3 D + (36^ + 234z/ 4 z/ 2 — 146z/2 

+ 261i/ 2 z/ 2 - 488z/ii/ 2 + 155 z/ 2 + 63- 155 v{ + 137zq - 585)}, (34) 

M i2 = —2 77 — — e(27ul + 135z/iz/( — 21v% + 135 z/V 2 — 210^i ^ 2 + 24z/ 2 + 27^ 

- 21/y 2 + 24^i - 155)77, (35) 

M13 = — -z/3 — —erq |4 a/3(9z/iZ / 2 + 10z/ 2 + 9uf — 6z/ 4 — 4)77 — (216z/ 2 + 288z/ 4 z/ 2 

- 154z/ 2 + 171i/ 2 - 38z/, + 420)}, (36) 

M14 =-—— 1/3 + — £i'' 3 }4(18z/2 + 27^i z/ 2 — 2z / 2 + + 14iq — 12)77 + \/3(51r / 2 

+ 114zqi/ 2 + 2 i / 2 + 87z/ 2 - 120i/i + 155)}, (37) 

M 2 1 = 277 — —e |4(9z/f + 45 z/iz / 2 + 9i/| + 45z/z/ 2 — 30uiu 2 — 18i/ 2 + 9z/f + 9i/ 2 

- I81/1 + 61)77 + 3a/3 1/3(12^2 - 6zqz/ 2 + 14i/ 2 - 15i/ 2 + 4zq - 5)}, (38) 

M 22 = 77~ — —^(Sz^ + 12uiu 2 — 18z/ 2 + 3^1 — 18^1 + 10)77 2 — 3\ // 3[i , i — ^ 2)^3 

x (9 z / 2 + 9z/i+4)77}, (39) 

M 2 3 =-——^3 + —£i^3}4(18z/2 + 27zqz/ 2 + 8 z / 2 + 9^ + 16zq — 12)77 + a/3(36 z / 2 

+ 72z/iz/ 2 - 54z / 2 + 81^ 2 - 90zq + 160)}, (40) 
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M 24 = ~z/ 3 + -^£i/3|4-\/3(9z/iz/2 + 8u 2 + 9z/ 2 - 4 )D - 9(21z/ 2 + 14 //i// 2 - 10z/ 2 

+ 13^—8//, +45)}, (41) 

M 3 1 = (Z) 2 - 3) - is|4V3 (i/ 3 - i/i)(7 - 9// 2 )// 2 Z) + (36z/f + 234^^ - 146z/ 3 2 

+ 261z/ 2 z/ 3 - 488z/!Z/ 3 + 155z/ 3 + 63//? - 155z/ 2 + 137/q - 585)}, (42) 

M 32 = —27)-e(27z/f + 135z/iz/| — 21z/| + 135z/ 2 z/ 3 — 210 z/i// 3 + 24z/ 3 + 27/zj 5 


21 z/ 2 + 24i/]. - 155)7), 


(43) 


M 33 = Z) 2 


- 3 + ~z^2 - |4\/3z/ 2 (9z/f + 9zqz/ 3 + 8z/ 3 - 4zq - 4)7) - (180z/f + 270zqz/| 


- 224z/ 2 


+ 198z/ 2 z/ 3 + 8 z/iz/ 3 + 419z/ 3 + 108z/f — 54z/ 2 + 3217^1 + 165)}, 


■W34 — I 2D + 


3 ^ 


z/ 2 J + ^t£ |4(27z/f - 39/y. 2 - 27ufu 3 + 165z/ii/ 3 - 54z/ 3 + 36z/ 2 

- 102i/i + 191)7) + 3v / 3z/ 2 (51z/| + 114^z/ 3 + 2z/ 3 + 87z/ 2 - 120z/i + 155)}, 
z/| + 45z/iz/| + 9z/f + 45z/ 2 z/ 3 — 30z/i// 3 — 18z/ 3 + 9 v\ + 9z/ 2 


M4 1 — 27) — —£ 


18z/i + 61)7) — 3 a/3z^2(1 2z^ 3 — 6//]z/ 3 + 14z/ 3 — 15z/ 2 + 4zq — 5) >, 


M 42 = D 2 — —£ j(3z/f + 12z/ii/ 3 — 18 z/ 3 + 3z/ 2 — 18z/i + 10)7) 2 — 3Vs(u 3 — zq) 


x (13 — 9z/ 2 


M43 = 2D - —1/ 2 + — e 
4 32 


(44) 


(45) 


(46) 


(47) 


z/f — 19z/| — 9z/ 2 z/ 3 + 27v\v 3 — 2v 3 — 2v\ — \i)u\ — 49)7) 


+ v / 3(72z/| + 54 z/!Z/ 3 - 12z/ 3 + 36//( - 78//! + 145 )i/ 2 }, (48) 

9 1 r /- 

M 44 = 7) 2 -z/ 2 -£< 4(3// 2 + 12 z/i// 3 — 18// 3 + 3z/ 2 — 18z/i + 10)Z) 2 — 12v3(9z/| + 9 //i// 3 

4 96 L 

+ 12// 3 - 4//! - 4)i/ 2 Z) - 27(21 v\ + - 10z/ 3 + 13z/( - 8u 3 + 45)z/ 2 }. (49) 

In a similar manner to the Newtonian case, we obtain the secular equation at the 1PN 
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order as 


A 2 A 6 + 2 <j 1 - A £ (77 _ i 0 y) }> A 4 + <j 1 + - ^£(308 + 12651/ + 162IF - 378 V 2 ) \ A 2 

+ T { v ~ 54 £ ( 521V ~ 72W ~ 126V ^ 


where W = viv 2 v^. 

Neglecting the trivial root A = 0, we obtain a cubic equation of r = A 2 as 


= 0, 
(50) 


t 3 + ar 2 + (3t + 7 = 0, 


(51) 


where 


a = 2 jl-ie(77-10y)j, 

P = 1 + (308 + 12651/ + 162VF - 378V" 2 ) , 


7 


1 16 

27 f 1 

V - —£(5211/ - 72 W - 126V 2 ] 


(52) 

(53) 

(54) 


4 ( 24 

In a similar manner to the Newtonian case, the PN triangular solution is stable if and only 
if all the roots r of Eq. (T5T1) are negative real, so that all A = ±\/t have no real part. 

In the 1PN approximation, the PN corrections to the Newtonian roots must be small. 
Thus, we can factor Eq. (T5TT) as 


(r + 1 — ae){r 2 + br + c) = 0, 


(55) 


where a, b, and c are constants and the 2PN terms are neglected. From Eqs. (I5TD and (1551) . 
one can obtain 


a = \-A77V - 14V 2 - 36 W), 

8 v 

(56) 

6 = 1 — -^-(771/ - 6 V 2 + 36VF)£, 

8 v 

(57) 

c = —V - — (1305F - 378V 2 + 162 W)e. 

4 16 

(58) 


Since £ C 1, we have — 1 + ae < 0, b > 0, and c > 0. The roots of Eq. (155]) are expressed as 

-b ± Vb 2 - 4c 


Ti — —1 + ae , r± = 


(59) 
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Hence, Eq. ()55]) has three negative real roots if and only if all the roots (1591) are negative 
real. Namely, it is necessary and sufficient for the stability that 

b 2 -AO 0. (60) 


For the critical value b 2 — 4c = 0, V in the PN terms of Eq. (1601) with Eqs. (157)1 and 
(l58]l can be replaced by the Newtonian critical value as V = 1/27, because 1PN corrections 
to it make 2PN (or higher-order) contributions and they can be neglected. Therefore, the 
condition for stability of the PN triangular solution becomes 




This is explicitly rewritten as 

miTTT.2 + 777.2777.3 + r773?77i 


+ 


15 


777x77727773 


1 ft 391 
11 ~ l4 e 


(61) 


(62) 


(t77i + 7772 + 777 3 ) 2 2 (t77i + 777 2 + 7773) 3 27 

Equation (162]) recovers the Newtonian condition ([[]) in the limit e —> 0. The PN correction 
in the right-hand side of Eq. (I62j) is negative and the PN term of the triple product of 
masses in the left-hand side of Eq. (162]) is positive. Hence, the PN condition for stability 
is tighter than the Newtonian one for any positive small value of the parameter e. Figure 
[3] shows the Newtonian stability regions Eq. ([TJ) and the 1PN ones Eq. (j62j) when e = 0.01 
(i.e. the order of magnitude of the PN effects is 0.01), for instance. 

Finally, we focus on the restricted three-body limit as n 3 —* 0 (i.e. W —> 0). In this case, 
the stability condition Eq. (I60|) becomes 

77- 1311U + 378U 2 


1 — 27V 


This is a quadratic inequality of V. Solving Eq. 


777-1777-2 J_ . 

(777! + m 2) 2 < 27 


-£ > 0 . 

for V, we obtain 
391 


(63) 


(64) 


where we used a relation 0 < V < 1/4 in the restricted three-body problem. Using a relation 
v\ + U 2 = 1 and assuming v\ > V 2 without loss of generality, we can rewrite Eq. (|64li as 

7772 17\/69 . . 

—;— < ^0 —( 65 ) 
777 1 + 7772 486 

re Newtonian value /i 0 = (9 — \/69)/18. This condition is in agreement with previous 


where t 
results 


24 


26] (see Eq. d3j) in the present paper). Figureshows a region of stability in the 


PN restricted three-body problem. 
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IV. CONCLUSION 


We examined the PN effects on the stability of the triangular solution in the relativistic 
three-body problem for general masses. The stability regions still exist even at the 1PN 
order. The PN stability condition Eq. flb2ll is consistent with the Newtonian one Eq. (JT]) 
in the limit as £ —> 0. The PN correction in the right-hand side of Eq. (|62|) . which is in 


agreement with that in the PN restricted three-body problem 


24 


26], makes the condition 


more strict than the Newtonian case for any small positive value of the parameter e. 

The PN term of the triple product of masses in the left-hand side of Eq. (l62jl appears 
not in the restricted case but in the general one. The instability is also enhanced by this 
term, while this effect is smaller than the other PN one in the case of mass ratios for stable 
configurations. If a system is mildly relativistic as £ = 0.01, for instance, the maximum 
value of W is O(10 -4 ) when v 2 = z/ 3 ps 0.019 in a stability region. Namely, the contribution 
from W is comparable to the 2PN (or more higher) order. This implies that triple systems 
with the PN triangular configuration for three finite masses are possible as well as restricted 
three-body systems. 

The PN triangular configuration ought to emit gravitational waves 7, [10]. Such a system 
will shrink by gravitational radiation reaction if its configuration is initially stable, and the 
PN effects on the long-term stability should be incorporated. In addition, we concentrate 
on the circular orbit in this paper. Gravitational radiation is known to decrease eccentricity 
in binary orbits. It is left as a future work to study gravitational radiation reaction to the 
PN triangular configuration. 
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Appendix A: A Derivation of the Equations of Motion for Perturbations 

We consider four perturbations in the orbital plane (see Fig. ED. First, we put the 
distances between the bodies as 


ru = £(1 + pu + Sxij), (Al) 

where Xij(— Xji ) is a perturbation in the distance 77 j and S is a bookkeeping parameter 
that denotes the smallness of the perturbations. By these perturbations, each angle <pu 
between vki and r jk (/ ^ J ^ K ) of the PN triangle is changed as 

1 Pij = — + V3p u + (A2) 

The perturbations Xij and Vhj relate to each other through the cosine formula, and the 
number of independent perturbations is three. 

The remaining one of the degrees of freedom corresponds to a change in the angular 
velocity of the bodies: 

6ij = ©/j + (A3) 


where 9jj and 07 j denote the direction of rjj to the reference frame and a perturbation in 
it, respectively. Qjj is the unperturbed direction which satisfies the equation as 

dQu 


Differentiating relations as 


dt 


— (Un(1 + wpn)- 


(A4) 


we obtain 


$23 — $12 — 7T ~ <^31; 
$31 = $12 + 7T + <^23, 


(AS) 

(A6) 


Dcr 23 = D(cr - -0 31 ), 
Da 31 = D(a + 1/7,3), 


(AT) 

(AS) 


where D denotes a differential operator with respect to a normalized time t = uj-^t and we 
denote a ±2 simply as a. Thus, the number of degrees of freedom for ( 07 2 , 073 , 031 ) is one and 
it corresponds to a change in the angular velocity. 
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Note that the perturbations have not only the Newtonian terms but also the 1PN ones. 
For instance, the perturbation er can be expanded as 


O' — (Tn + Upn, 


(A9) 


where <tn and opn(= 0(e)) are the Newtonian term and the 1PN one, respectively. In the 
following, we neglect the terms of second (and higher) order in 5. Namely, we calculate to 
the terms of order e x 5 (i.e. the linear perturbation at the 1PN order). 

Using a complex plane as the orbital one, we denote the relative position of the bodies 
as rij —> Zjj = rue ldlJ . The EIH equation of motion for z\ 2 becomes 

d 2 z V2 


dt 2 


= F 12 e 


iOi2 


(A10) 


The left-hand side of this equation is 
d 2 z\ 2 

= [~~{1 + 2o)pn + Pv 2 + d(2Da 12 + X 12 — D 2 \ 12 + 2u)pnF><7 12 + 2p 12 Dcr 12 

+2cupnXi2)} + id(2Dx\2 + F) 2 a\ 2 + 2uj p ^Dx\2 + Pi2-D 2 <7i 2 )] e ldl2 . (All) 

F\ 2 in the right-hand side of the equation of motion can be expanded as 


F 12 — Tn 12 + eF PN12 + 5pNper 12 + £^F] 


PNper 12) 


(A12) 


where Fni 2 and Fpn 12 are the unperturbed Newtonian and PN terms, respectively, and 
FN P eri 2 and FpNper 12 are the perturbed Newtonian and PN terms, respectively. These are 

M 

Fn 12 = — -p, (A13) 

Fpn 12 = ^^-(45z /1 + 54zqz/ 2 - 60z/ 2 + 45z/( - 60z/i + 97), (A14) 

1 M. , , . , ,3a/3 M . , , A . 

FN P eri 2 = ^~jp [3^3(A23 + X 31 ) + 2(2 — 3z/ 3 )xi 2 ] + -^-^{xsi — X 23 ), (A15) 
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F, 


PNper 12 


1 M 
161 ? 


—2(54^2 + IO 8 Z/ 1 Z /2 — 86^2 + 108^Z2 2 — 82z/iz/ 2 + 167z/ 2 + 54^ 


— 86 ^ + 167zq — 29)xi2 — ^ 3 ( 45^2 + 108z/iz/ 2 + 8 z / 2 + 90z/^ — 108z/i + 150) 
X X 23 - Us(90^2 + 108iqz/ 2 - 108z/ 2 + 45^ + 8 zq + 150)y 3 i + 8(^2 - 

+ 2z^2 — ^ 1^2 — 4zqz/ 2 — 7z/ 2 + + 2z/^ — 7ui)Dai2 — 2v 3 (y2 4- 22zqi/' 2 + 4z/ 2 

+ 4z/^ + 2zq + 8 )-Dct 23 — 2z/ 3 (4z/ 2 + 22/qiq + 2z/ 2 + + 4iq + 8 )Z)c 3 i 

+ 8\/3i/ 3 (z2i - i/ 2 )(3 - z/ 3 )Fxi 2 - 2a/3z/ 3 (9z/ 2 - 4zqz/ 2 - 4z/ 2 + 4u\ + 6zq 

- 16)F>x 23 + 2 v / 3 z/ 3 ( 4 z /2 - 4zqz / 2 + 6 u 2 + 9zq 2 - 4zq - 16)Fx3i 


+« 


1 M 
16~^ 


-12v / 3i/ 3 (z/i - z/ 2 )(3 - z/ 3 )xi 2 + V3v 3 (57v% + 36zqz/ 2 - 24z/ 2 


+ 42 — 12ui + 130)y 23 — \/8v 3 {42v\ + 36zqz/ 2 — 12z/ 2 + 57^ — 24iq 

+ 130)x 3 i — 8\/3z2 3 (^i — UjXul + V 2 )D(Ji 2 + 2\^3u 3 ( k i i 2 — 12iqz/ 2 + 14z/ 2 

— 4z^ + 10zq + 8 )D<t 23 + 2a/3i/ 3 (4z/ 2 + 12 zqz/ 2 — 1(V 2 — — 14 zq — 8 )F<t 3 i 

— 8(3z/2 + 9z/iz/2 — 6z/2 + 9^i^ 2 — 8zqz/ 2 — 5z/ 2 + 3/^f — 6//f — 5zq)77;\q2 

— 2z/ 3 (9z/2 + 30zqz/ 2 + 10z/ 2 + 12^ — 18zq — 16)77x23 — 2v 3 (12i>1 + 30zqz/ 2 


- 18zq + 9^i + 10zq - 16)Z7 x 3 i 


(A16) 


For exchanging indices between 1 and 2 such as zq yy z/ 2 and X 31 -B- X 23 , we have the symme- 
try/antisymmetry in the real/imaginary parts of Eqs. (1A13j) - (1A16D . This (anti)symmetry 
may ensure the form of the equations, though they are rather complicated. We can obtain 
the EIH equations of motion for z 23 and z 3 i by the cyclic manipulations asl —> 2 —)■ 3 —> 1. 
Since the unperturbed terms in the equations give the PN triangular equilibrium solution, 


we focus on the perturbed terms. It is convenient to transform the variables as 

1 


23] 


X23 — 2 


(1 - 3pi 2 )x 3 i + (1 - 3p 3 i)xi2 + v^3(l - p23)^23 

X = X 31 — Xl 2 - 


(A17) 

(A18) 


Using the relations Eqs. (IA7j) , (1A8I1 . (1A17|) . and (IA18p . we obtain the equations of motion 
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m 3 



FIG. 1: PN triangular configuration. Each body is located at one of the apexes, pu denotes the 
PN corrections to each side length at the 1PN order. In the equilateral case, pi 2 = P23 = P31 = 0, 
namely, ri 2 = 7-23 = ^31 = £ according to Eq. GDI- 
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FIG. 2: Four perturbations in the PN triangular configuration. The perturbations X12 and X 31 
denote changes in r \2 and 7 - 31 , respectively, ^23 is a perturbation in the interior angle ^ 23 , and a 
corresponds to a change in the angular velocity. 
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FIG. 3: Mass ratios for a stable configuration that satisfies the Newtonian condition Eq. (jT]) and 
the PN one Eq. (1621) when e = 0.01 (i.e. the order of magnitude of the 1PN effects is 0.01), for 
instance. For values of the mass ratios within the colored areas, the triangular configuration for 
three finite masses is stable. Top: All the stability regions. Bottom: The regions around small 
u\ and i/ 2 , where the third mass is dominant. The stability regions at the 1PN order still exist, 
though they are more narrow than the Newtonian case. 
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FIG. 4: The stability region Eq. (1651) for the PN restricted three-body systems in the (y 2 , e) plane. 
For values of V 2 and e within the colored area, the PN triangular configuration in the restricted 
three-body case is stable. This figure corresponds to Fig. 1 in Ref. 24]. 
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